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Abstract 

Classical and nonclassical symmetries of the nonlinear heat equation 

Ut^Uxx+f{u), (1) 

are considered. The method of differential Grobner bases is used both to find the conditions 
on f{u) under which symmetries other than the trivial spatial and temporal translational sym- 
metries exist, and to solve the determining equations for the infinitesimals. A catalogue of 
symmetry reductions is given including some new reductions for the linear heat equation and a 
catalogue of exact solutions of (1) for cubic f{u) in terms of the roots of f{u) = 0. 
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1 Introduction 



The nonhnear heat equation 



Ut = Uxx + f{u), 



[l.D 



where x and t are the independent variables, f{u) is an arbitrary sufficiently differentiable function 
and subscripts denote partial derivatives, arises in several important physical applications including 




), in the 
ill) 



microwave heating (where f{u) is the rate of absorption of microwave energy, cf. [pO| , 
theory of chemical reactions (where f{u) is the temperature dependent reaction rate, cf. 
and in mathematical biology (where f{u) represents the reaction kinetics in a diffusion process, cf. 

The classical method for finding symmetry reductions of pdes is the Lie group method of infin- 
itesimal transformations |Tl|, 49, 73, 78, , 93|. To apply the classical method to the nonlinear 
heat equation (IT) we consider the one-parameter Lie group of infinitesimal transformations in 
(x, t, u) given by 



X = X + £^{x,t,u) + 0{e^)^ 
t = t + eT{x,t,u) + Oie"^), 
u = u + £(t){x^t,u) + 0{e'^)^ 



[1.2) 



where e is the group parameter. Requiring that (1.1) is invariant under this transformation yields an 
overdetermined, linear system of equations for the infinitesimals ^(x,t, m), r(x,t, m) and 0(x,t,u). 
The associated Lie algebra of infinitesimal symmetries is the set of vector fields of the form 



^, . d , . d 
C(a;, t,u)— + t{x, t,u)— + 



(1.3) 



Classical symmetries of the nonlinear heat equation (1.1) have studied by Dorodnitsyn |2j] who 
classified conditions on f{u) for which symmetries exist, though he did not give the associated 
reductions (see also [^, ^). 

There have been several generalizations of the classical Lie group method for symmetry reduc- 
tions. Ovsiannikov [ |78| developed the method of partially invariant solutions. Bluman and Cole 
1^, in their study of symmetry reductions of the linear heat equation, proposed the so-called non- 
classical method of group-invariant solutions (in the sequel referred to as the nonclassical method), 
which is also known as the "method of conditional symmetries" |33, 34, 36, 37, 59, 81, |8^, |103f| 
and the "method of partial symmetries of the first type" | |97[ |. In this method, the original pde 
(1.1) is augmented with the invariant surface condition 



^p = S,{x, t, u)ux + t{x, t, u)ut — (j){x, t, u) = 0, 



(1.4) 



which is associated with the vector field (1.3). By requiring that both (1.1) and (1.4) are invariant 
under the transformation (1.2) one obtains an overdetermined, nonlinear system of equations for 
the infinitesimals ^{x,t,u), T(x,t,u) and (j){x,t,u), which appear in both the transformation (1.2) 
and the supplementary condition (1.4). Since the number of determining equations arising in the 
nonclassical method is smaller than for the classical method and since all solutions of the classical 
determining equations necessarily satisfy the nonclassical determining equations, the solution set 
may be larger in the nonclassical case. For some equations, such as the Korteweg-de Vries equa- 
tion, the infinitesimals arising from the classical and nonclassical methods coincide. It should be 
emphasized that the associated nonclassical vector fields do not form a vector space, still less a Lie 
algebra, since the invariant surface condition (1.4) depends upon the particular reduction. Subse- 
quently, these methods were further generalized by Olver and Rosenau |^, 76 1 to include "weak 
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symmetries" and, even more generally, "side conditions" or "differential constraints", and they 
concluded that "the unifying theme behind finding special solutions of pdes is not, as is commonly 
supposed, group theory, but rather the more analytic subject of overdetermined systems of pdes" . 

In this paper we consider both classical and nonclassical symmetries of the nonlinear heat 
equation (1.1), and find conditions on f{u) under which symmetries, other than trivial translational 
ones, may exist. The method used to find the necessary conditions on f{u), and from there to find 
the symmetries, is that of differential Grobner bases (dgbs) (6^, This method yields the 

"triangulation" of a nonlinear system of PDES from which solutions can be obtained in a systematic 
way. 

The basic algorithm used to generate a dgb for a system of pdes is the Kolchin-Ritt algorithm, 
which forms an essential part of the StandardForm algorithm 84 1 and the diffGbasis algorithm 



1 63]. It appears to have been formally written down first by Carra-Ferro [1^. The Kolchin-Ritt 
algorithm has been used for some time, intuitively if not formally, to yield integr ability conditions 
for linear and orthonomic systems With hindsight, one can see that the integr ability 

conditions generated by the Janet-Riquier algorithm, which completes a system of orthonomic 



pdes to involutive form (cf. for example |52, 53, 94, 95| and references therein), are produced 



by exactly the calculations performed in the Kolchin-Ritt algorithm. (It should be noted that 
several authors incorrectly assume that when the Kolchin-Ritt algorithm terminates, the output 
system is involutive, or passive. A counter-example can be found in [^.) The production of 
integrability conditions requires a total degree ordering on the derivative terms. On the other 
hand, the elimination ideals, or triangulation, of a system of pdes requires a purely lexicographical 
ordering. 

Recently, the limits of the Kolchin-Ritt algorithm in calculating a dgb of a nonlinear system were 
found, and sufficient conditions for obtaining a DGB for a nonlinear system of pdes were formulated 
1 61]. The Kolchin-Ritt algorithm generates a dgb only "up to" (in a sense to be defined) a certain 
set of differential coefficients, and these coefficients must not lie in the ideal. Proving that these 
coefficients do not lie in the ideal has no algorithmic solution as yet, and there are reasons as to 
why it may not have one at all. The condition is also one of the sufficient conditions in ]^^ for a 
DGB of a nonlinear system. Thus far, attempts to formulate algorithms that circumnavigate the 
difficulty, for example by calculating the algebraic Grobner basis of successive prolongations of the 
system |l^, lead to processes that do not terminate in finite time. Fortunately in many 

cases one can prove the condition directly, and there is a range of techniques that can be used when 
the condition fails. Despite the limitation, the Kolchin-Ritt algorithm is the fundamental tool in 
the theory, and improvements to efficiency and practicality of this algorithm is the most important 
problem to be faced from an applied point of view. A promising approach to this problem can be 
found in the thesis of Lisle J^] , who rewrites the systems of PDES using first order noncommutative 
operators to effect an improvement in efficiency. 

The next section contains an outline of the theory of differential Grobner bases. After establish- 
ing the notation and formulae required, we show how to calculate the triangulation of a system of 
PDES using this theory. A variation of the Kolchin-Ritt algorithm is given that appears to be much 
faster than the one presented in 63], at least for the examples calculated in this paper, and 
has been developed especially for use in solving the determining equations arising from the non- 
classical method. Such systems of determining equations appear to have a natural ordering which 
can be exploited to improve the efficiency of the calcualtions. We also give a "skeletal" form of the 
Kolchin-Ritt algorithm which we found useful in extracting information to reduce the complexity 
of the calculations. We refer to this second algorithm as "Direct Search" , which appears to be new. 

In the third section of this paper, classical symmetries of the general nonlinear heat equation 
are catalogued for all those nonlinear f{u) where a non-trivial symmetry is possible. Similarly, 
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nonclassical symmetries are given in the fourth section, including some nonclassical symmetries for 
the linear heat equation. A catalogue of nonclassical symmetry reduction solutions, for cubic /, is 
given in the fifth section; the classification is in terms of the roots of f{u). The paper concludes 
with some plots of several exact solutions of (1.1) derived from nonclassical symmetry reductions. 

There is much current interest in the determination of symmetry reductions of PDES which 
reduce the equations to odes. One frequently then checks if the resulting ode is of Painleve type 
(i.e., its solutions have no movable singularities other than poles). It appears to be the case that 
whenever the ODE is of Painleve type, or can be transformed to one that is, then it can be solved 
explicitly, leading to exact solutions to the original equation. Conversely, if the resulting ode is 
not of Painleve type, then often one is unable to solve it explicitly. 

Since solutions of pdes asymptotically tend to solutions of lower-dimensional equations ob- 
tained by symmetry reduction, some of these special solutions will illustrate important physical 
phenomena. In particular, for reaction-diffusion equations such as (1.1), exact solutions arising from 



symmetry methods can often be effectively used to study properties such as "blow-up" |^, g^, |46|| . 
Furthermore, explicit solutions (such as those found by symmetry methods) can play an impor- 
tant role in the design and testing of numerical integrators; these solutions provide an important 
practical check on the accuracy and reliability of such integrators (cf. |Q, |8^, |9H). 

Several exact solutions of equations of the form ( |l.lD have been derived in the literature. 
Ablowitz and Zeppetella [E| obtained an exact travelling wave solution of Fisher's equation 



Ut 



+ u(l 



(1.5) 



by finding the special wave speed for which the resulting ode is of Painleve-type. Recently Guo 
and Chen [41| have used the Painleve expansion method |100| , 101 ] to obtain some heteroclinic 
and homoclinic solutions of ( |1.5D . Kaliappan |54| and Herrera, Minzoni and Ondarza |48| have 
derived travelling wave solutions of the generalized Fisher's equations 



Ut 
Ut 



+ u- 



2p~l 



(1.6) 
(1.7) 



respectively, 
equation 



Cariello and Tabor |jT^, 13 1 found an exact solution of the real Newell- Whitehead 
(or Kolmogoroff-Petrovsky-Piscounov equation fS^) equation 



Ut 



Ux 



+ u{l 



(1.8) 



using a truncated Painleve expansion and verified that it derives from a nonclassical symmetry 
reduction (see also [^]). Several authors have studied exact solutions of the Fitzhugh-Nagumo 
equation 

Ut = Uxx + u{l - u){u - a) (1.9) 

where a is an arbitrary parameter, which arises in population genetics ^ and models the trans- 
mission of nerve impulses [ 30 , ^] . Travelling wave solutions of the Fitzhugh-Nagumo equation ( |1.9D 
have been studied by several authors 0, |2^, |9^ . Exact solutions of (|l.9D have been obtained 
using various techniques including Vorob'ev (who calls the associated symmetry a "partial 
symmetry of the first type"), by Kawahara and Tanaka using Hirota's bi-linear method j50[ |, 
by Hereman [45| using the truncated Painleve expansion method (see also |l^) and more recently 
by Nucci and Clarkson [^] using the nonclassical method (see also [^l|). Exact solutions of the 
Huxley equation 

Ut = Uxx+U^{l-u) (1.10) 
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have been obtained by Chen and Guo , using a truncated Painleve expansion (see also G% ) , 



and by Clarkson and Mansfield |21| using the nonclassical method. Exact solutions of equation 

(1.11) 



ut 



have been obtained by Clarkson and Mansfield [^] using the nonclassical method. 

A number of symbolic manipulation programmes have been developed, e.g. in macsyma |16], 

8|, |3, |5l, MATHEMATICA [|, 0, MUMATH Q and REDUCE jS^, |9|, IS^, ||, |9|, IDH] , that 



MAPLE 



calculate the determining equations for classical Lie symmetries of pdes. A survey of the different 
packages presently available and a discussion of their strengths and applications is given in |46]. 
The determining equations for symmetries solved in this paper were calculated using the macsyma 



programme SYMMGRP.MAX 
the MAPLE package diffgrob2 



The triangulations of these systems of pdes were calculated using 
62|. The plots of some of the exact solutions of (1.1) derived here 



were drawn using maple and mathematica. 



2 Differential Grobner Bases 

In this section we discuss the triangulation of a system of pdes from which solutions may be obtained 
in a systematic way and give some algorithms using which a triangulation can be determined. First 
we introduce some required notation. 

2.1 Some Notation 

The determining equations for symmetries of PDES can be regarded as polynomials in the infinitesi- 
mals ^, T and (p, their derivatives, and the variables {x, t, u}, with complex coefficients. We establish 
our notation and formulae required for the statement of the algorithms in generic co-ordinates for 
ease of reference. 

Let / be a differential polynomial; that is, a polynomial with complex coefficients in the inde- 
pendent variables {xi, . . . , the functions {ui, . . . , Um} and the derivative terms {D°'uj} where 

D°u- - ^ 

'~ . . . dx'^" ■ 

The set of all such differential polynomials (d.p.'s) is denoted Rn,m or 

, Xji , Ul ) • • • ) "Wm ) ■ 

All the concepts depend upon an ordering on the derivative terms. There are many suitable 
orderings for different purposes. They need to satisfy B"uj > D^iifc implies that D'^D^u^ > D^D^Ufe 
and ^^Uj < T)'^T)°'uj. We assume that Uj > x^ for all j and k. The standard lexicographic ordering 
based on ui < U2 < • • • < Um and xi < X2 < . . . < x„ is given by 

if Uj > Uk , 

else j = kon > l3n, 

else On = f^n, ■ ■ ■ , dn-j = /?n-j; On-j-i > Pn-j-i for some j such that < j < n — 2. 

The highest derivative term occurring in a d.p., /, is denoted HDT(/). 

The highest power of the HDT(/) occurring in / is denoted Hp{f). 

The highest coefficient, HcoeS{f), is defined to be coeff(/, HDT(/)^*'(''^)). 
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The head of / is Head(/) = i/coeff(/)HDT(/)^P(^). 

The separant of /, the highest coefficient of D"/ for any non-zero multi-index a, is denoted 
Sepif). 

A pseudo-reduction of a d.p. / by G = {gi, g2, ■ ■ ■ ;5fc} C Rn,m effects ehmination from / of any 
terms that are of the form D"HDT(5(j), for some i G {1, . . . , fc} and a G N"". 

Let a derivative term DT occur in / to some power p. Suppose there exists an a G such that 
D°HDT((7) = DT for some g & G. If a = assume further that p > Hp{g). A pseudo-reduction to 
/' of / by 5 is given by the formulae 



gcoeff(D"ff)/ - coeff(/, DTP)DTP-'B-g ^ H a ^ 0, 
ffcoeff(g)/ - coeff(/, £)rP)HDT(g)(^-^^(g))g ^ ^ ^ 



where Z is the greatest common divisor of HcoeS{D°'g) and coeff(/, DTP). 

The pseudo-normal form of / with respect to a set G is obtained when no further pseudo- 
reduction with respect to any member of G is possible, and is denoted normaF(/, G). 

The diffS polynomial of two d.p.'s is obtained first by cross-differentiating and then by cross- 
multiplying by the highest coefficients, and subtracting. 

Let /i and /2 be two d.p.'s with the same highest unknown. Let ai and 0:2 be the smallest 
multi-indices possible such that D"iHDT(/i) = D"2HDT(/2). 

Let Z = gcd(i/coeff(D"Vi), i/coeff(D"2/2)). If both ai, 02 7^ 0, define 

w, -•^ gcoeff(D"Vi)D"V2 - gcoeff(D"V2)D"Vi 
diffS'poly(/i,/2) = . 



If ai = and 02 / 0, then 

^^ gcoeff(/i)HDT(/i)(^P(/^)-^)D"V2 - gcoeff(D"V2)/i 
ditt6poly(/i,/2) = ^ 

and similarly if cti / and 02 = 0. 

If Q;i = 0:2 = so that HDT(/i) = HDT(/2), or if /i and /2 have different highest unknowns, 
then the differential S polynomial is defined to be 

^■ffc ^ ^f f\ Head(/2)/i - Head(/i)/2 

dlff5poly(/i, f2) = ,/TT n/ . N TT ^ff ^^ • 

gcd(Head(/i), Head(/2)) 

For linear equations with different highest unknowns, we take their diffSpolynomial to be zero. 

For G C i? n,m let S{G) be the multiplicative set in Rn,m generated by the set of factors of all 
the highest coefficients and separants of the g in G. We assume that {1, —1} G S{G). 



2.2 The Triangulation of a System of PDEs 

Given a system of PDES, S C Rn,mt we define the ideal generated by S to be all those equations 
that can be obtained from the elements of S by differentiating and adding, and multiplying by 
arbitrary elements of Rn,m- In other words, 

/(S) = {5^ga,iD°/, : G S, ga,^ e Rn,m, oc G N"}. 

If the system is linear, we allow multiplication by polynomials in the variables xi, . . . ,x„ over C 
only. 
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A differential Grohner basis of /(S) is defined to be a set of generators G of /(S) sucli that 
every element of /(S) pseudo-reduces to zero with respect to G. 

Differential Grobner Bases (dgbs) of ideals of differential polynomials have many useful prop- 
erties and can be used to solve, in theory at least, a great many problems. These include being able 
to find all elimination ideals, integrability conditions and all compatibility conditions of a system 
of nonlinear PDES [63]. 



A DGB of a system of pdes depends on the ordering of derivative terms. Given a dgb for a 
system E of nonlinear PDES, we have the following theorem [6^ : 

THEOREM: Elimination Ideals. If S is a dgb for /(E) in a lexicographic ordering with ui < 
U2 < ■ ■ ■ < Urn and xi < X2 < ■ ■ ■ < Xn, then T,p^q := S n Cdifr(2;i, . . . , Xp, ui, . . . , Uq) generates, up 
to S'(S) Ip^q = /(S) nCdifj(xi, . . . , Xp, ui, . . . ,Uq). That is, for every f E Ip^q there exist an s £ S{T,) 
such that sf £ 'Sp,q that is, sf is in the ideal generated by ^p^q- ° 

In other words, every condition involving only the first q unknowns with derivatives with re- 
spect to the first p variables, obtainable from the original set of pdes by differentiating, adding, 
multiplying by like d.p.'s and so on, can be "read off' from a dgb of the system in a lexicographic 
order, up to S{T,), which represents the "margin of error" in the theory. 

For systems S that are linear in their highest derivative terms, the set S{T,) is trivial. The 
difficulties posed by non-trivial S{T,) are caused by the use of pseudo-reduction in which one is 
allowed to multiply the expression by non-trivial coefficients. However replacing pseudo-reduction 
by reduction, in which only multiplication by constants or polynomials in the independent variables 
is allowed, leads to an algorithm that will not terminate on general nonlinear systems. Even 
if terminated artificially, algorithms using reduction have two additionally problems. First, the 
expression swell renders the method impractical even for small problems, while second, there is 
currently no method for deciding at which level of prolongation of the system to terminate. 

The sets /(S) n Cdiff(xi, . . . , Xp, ui, . . . , Uq) are called the elimination ideals of /(S), while the 
dgb listed in decreasing order, with respect to the standard lexigraphic ordering described above, 
is called the "triangulation" of the system 

2.3 The Algorithms 

The basic algorithm used to calculate a DGB is the Kolchin-Ritt algorithm 61, 33, |8^ . 
In this algorithm, pairs of d.p.'s are cross-differentiated and cross- multiplied to make the leading 
derivative terms in each d.p. cancel. The result is then pseudo-reduced with respect to all other 
d.p.'s to obtain a new d.p. This is done systematically until the result of cross-differentiating all 
pairs of d.p.'s given or obtained pseudo-reduce to zero. 

In practice, there are three difficulties in calculating a dgb. The first is the problem of the 
build-up of differential coefficients obtained by repeated cross-multiplication and pseudo-reduction. 
If the coefficient of a highest derivative term of any d.p. contains a derivative term, then Kolchin- 



Ritt no longer necessarily suffices to calculate a dgb. Sufficient conditions may be found in |63], 
but there remains the difficulty that one of the sufficient conditions given in [^] is not algorithmic. 
This condition is that none of the highest coefficients or separants of the equations comprising 
the system lie in the ideal they generate, that is, S{Ti) n /(S) = (cf., §2.2). In addition, if 
this condition fails then pseudo-reduction can lead to spurious zeroes in the calculation of the 
algorithms, resulting in loss of information. It is well known that the complexity of the Buchberger 
algorithm for computing the Grobner basis of a polynomial ideal is doubly exponential [^]. Since 
the Kolchin-Ritt algorithm is modelled on the Buchberger algorithm, we have that the complexity 
of the Kolchin-Ritt algorithm is at least that of Buchberger's. However the main problem is that 
of "expression swell" in which the length of the expressions calculated exceed the memory limits of 
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the available computer. 

In spite of these problems, a dgb can often be obtained by a clever choice of the order in which 
pairs are cross-differentiated. We present here a version of the Kolchin-Ritt algorithm which for 
some systems is much faster than that given in |15, This version of the Kolchin-Ritt algorithm 
sorts the system after each new condition is found so that we are essentially finding all the lowest 
conditions (relative to the order) first. This seems to work well with those systems that are already 
partially sorted with respect to some lexicographical ordering on the dependent and independent 
variables, such as those in this paper. 

The output statement of the algorithm shows the limits of the Kolchin-Ritt algorithm when one 
begins with, or obtains en route, an equation where the coefficient of the highest derivative contains 
differential terms. The proof of the output statement can be found in |63]. It shows that when 
the highest coefficient of a d.p. vanishes, for example when considering additional constraints or 
conditions, then one must re-do the calculation taking the additional conditions into account from 
the beginning. In particular, when an equation obtained factorises, and one wishes to consider a 
factor that does not include the highest derivative term for that equation, then one must re-start 
the algorithm including the desired factor at the outset. 

We make some remarks about factors. It should be mentioned that if a product of functions is 
zero, say /1/2 = 0, then in general one can only conclude that one of /i, /2 must be zero if the fi 
are analytic. However, if one obtains, for example, a cubic in u^, then the continuity of Ux ensures 
that Ux must be one of the roots of the polynomial. 
Algorithm Kolchin-Ritt (sorting) 

INPUT: S a system of d.p.'s, a termordering 
OUTPUT: S' D S such that 

5(s') n /(s') = ^ V/ G /(s') 

3 s £ S{T,') such that sf pseudo-reduces to with respect to E'. 



pairsdone = {} 
i = l,j = 2 

sort S into increasing order 
S' = S 

while i / do 

let fk be the k^^ element of S' for < /c < 
if not {fi, fj} G pairsdone then 

pairsdone=pairsdone U {{fi, fj}} 

h =diSSpoly{fi, fj) 

h =normaF(/i, S') 

if /i / 0, then S' = S' U {h} 

sort S' into increasing order 
i = l,j = 2 
else if j = then j = i + 2, i = i + l 
else j = j + 1 

end 

We give now a related algorithm, the "direct search" algorithm. The "direct search" algorithm 
recursively looks for conditions that are lower in the order than those given. Its output is a "skele- 
ton" of that of Kolchin-Ritt, but often provides useful information that allows the triangulation of 
the system to proceed faster. 
Algorithm Direct Search 
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INPUT: F a set of d.p.'s, a termordering 
OUTPUT: Gd F 



G = F 

sort F into increasing order 

let fk be the k^^ element of F for < /c < |F| 

i = \F\-l 

h = f\F\ 
while i 7^ do 

while /i 7^ do 

G = GU {h} 

k = h 

h =normaF(diffSpoly(/i, /,), G) 
h = min{A;, /J 
i = i — 1 

end 

The direct search algorithm should be regarded as a "template" algorithm only. It is a natural 
generalization of the well-known process of setting coefficients to zero in the situation where one 
has an expression that is polynomial in {zz, Z2, Zm} say, with coefficients that are independent 
of the Zi. 



3 Classical Symmetries 

The determining equations for the classical symmetries of the nonlinear heat equation (1.1) are 



/l 


Tu = 0, 




/2 


Tx = 0, 




/3 


Tt - = 0, 




h 


C« = 0, 




h 


4>uu = 0, 




k 




= 0, 


fr 


4>t - (pxx + 4>uf 


- - Vix = 0. 



These equations were calculated using a macsyma programme symmgrp.max [16|. The spde 
package in reduce |87] appears unable to handle equations such as (1.1) which contain arbitrary 
functions. 

Assume the lexicographic ordering given hy (j)>T>^>f and x > t > u. From /2 and /s we 
immediately obtain 

/s : ixx = 0. 

Calculating diffS'poly(/6, /y) we obtain 

/g : -2(/)„t + 2(i)fuu + ^fuix - ixt = 0. 

Performing the "direct search" algorithm on /s and /g and reducing the results with respect to 
{iu,ixx,(l>uu} yields 



h : 
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^2 • ^ fuufuuuu ~l~ 2 (fuuu) ^ 4* ~^ "^fuufuuuix — 0, 

^3 • (,fuu) ^fuuuu (fuuu) '^fuu (fuuuu) ~\~ fuufuuufuuuuu^ — 0. 

Performing a "direct search" algorithm on /g and /g and reducing the resuhs with respect to 
{iu.ixxAuu} yields 



Ji 

32 



4^xufuu ~l~ 4xfuuu — 0) 
itt + '^(f>xfuu = 0, 
fuuuitt + (fuu)^ it = 0. 



Since = we must have, differentiating with respect to u and cross-multiplying, that 

tt ^(fuuu) fuufuuuu^ fuu 0- 

We now consider the possibilities for f(u); although it is possible that other "direct searches" 
will find other conditions, we can simplify the calculation of other conditions by putting = for 
example. 

We solve those factors in and j4 involving derivatives of f(u). Setting H = fuu in ^3 we 
obtain (assuming that H ^ and 7^ 0) 

Huu(Huf' — '2,H(Huu)^ + HHuHuuu = 0. 
Dividing by HHuHuu we have 

Hu 2 -^uu ^ Huuu g 

H Hu Huu 
This integrates easily; the possibilities from are 



(i) 


f(u) 


= av? + hv? + CU + d, 


(ii) 


f(n) 


= (au + fe)"""*"^ + CM + d, for n / 0, 1, 2, 3, 


(iii) 


f(n) 


= ae*'' + cu + d, 


(iv) 


f(n) 


= [ln(an + h)]/a + cu + d. 


(v) 


f(u) 


= [(au + b) ln(au + b) — (au + b)] / a^ + cu + d, or 


(vi) 


ix = 






where a, 6, c and d are arbitrary constants. 

Next we look at the equation j^. In this case we obtain 

(i) f(u) = au^ + bu + c, 

(ii) f(u) = [(au + b) ln(au + b) — (au + b)] / a^ + cu + d, or, 

(iii) it = 0. 

We shall not consider f(u) linear for the classical case as the symmetries in that case are already 
well known Q . Hence for those solutions f(u) not of the form u ln(n) + bu-\-c or au^ + 6u + c, with 
b and c constants, we must have that it = 0. 

If ix = 0, from the conditions ki and k2 one can conclude that if f(u) is not of the form 
It In (it) + 6tt + c or au^ + bu + c, then only the trivial solution, i = io, t = tq, 4> = with and tq 
constants, is possible. This solution exists for all choices of f(u) and corresponds to the travelling 
wave reduction 

u(x, t) = w(z), Z = X — ct 
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where c is an arbitrary constant and w{z) satisfies 

w" + cw' + f{w) = 0. 

The cases where f{u) is of the form uln{u) + bu + c or av? + bu + c require further attention 

since then the equation k2 does not provide a condition for (p. For the other possibihtics for f{u), 

we must re-calculate the triangulation to obtain further conditions on ^. For the remainder of this 
section, ki, K2, K3 and K4 are constants. 

Case 3.1: f{u) = + bu + c. Here the triangulation of the system is (in addition to /i, /2, /a, 
/4): 

4<P + 8uC^ + 46Cx + SCxt = 0, 6 = 0, ^xx = 0, Ub^-4c)=0. 

Thus unless f{u) is a perfect square, the only solution to the system is the trivial one. If f{u) = 
then 

^ = KiX + K2, T = 2Kit + K3, (f) = —2kiU. (3.1) 

Hence for 7^ (set ki = 1) we obtain the symmetry reduction 

U{x,t) = {t + ^K3)-^w{z), Z= {x + K2)/it + K3/2)^/^ (3.2) 

where w{z) satisfies 

w" + \zw' + w; + u;^ = 0. (3.3) 

It is easily shown that this equation is not of Painleve-type. 

Case 3.2: j(u) = uln(n) +6u-|-c. In the case c = 0, the triangulation of the system (in addition 

to /i, /2, /s, /4) is given by 

- ^ = 0, - ^ = 0, 2<\)^ + uit = 0, = 0, itt -6 = 

which yields 

^ = «;ie* + AC2, T = Ks, (f) = u(^—^Kix + Kij e^. (3.4) 
Hence for K3 ^ (set K3 = 1) we obtain the symmetry reduction 

u{x, t) = w{z) exp i^ — ^Kixe^ + ^K^e^^ + (^^K,iK2 + e^^ , z = x — Kie^ — H2t (3.5) 

where w{z) satisfies 

w" + n2w' + hw + w\nw = Q. (3-6) 

It is easily shown that this equation is not of Painleve-type. 

If K3 = (set Ki = 1 and K2 = 0) , we obtain the symmetry reduction 

„(x.«) = .(.)exp{ -^;f^;.r } , = e' (3.7, 

where v{s) satisfies 

s^^=v\n{v) + {h + nl-\)v (3.8) 

which has general solution 

v{s) = exp |/cs — 6 — + || (3-9) 
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with k an arbitrary constant. Hence we obtain the exact solution 

u{x,t) = expjfce* - \{x - 2K4f - b + ^} . (3.10) 

The solutions for the other cases of f{u) are calculated the same way and are listed below. 
Case 3.3: /(n) = n'^ + au^ + bu + c. Unless ab — 9c = and — 36 = there is only the 
trivial solution. The two algebraic conditions together imply that f{u) is a perfect cube. Setting 
f{u) = yields 

^ = KlX + K2, T = 2Klt + K2, (f) = —KlU. (^-ll) 

Hence for ki ^ (set aci = 1) we obtain the symmetry reduction 

u{x, t) = {t + K^/2)-^l'^w{z), z = {x + K2)/{t + n^/2f/'^ (3.12) 
where w{z) satisfies 

w" + \zw' + \w + w^ = Q. (3.13) 
It is easily shown that this equation is not of Painleve-type. 

Case 3.4: f{u) = + 6n + c, n 7^ 0, 1, 2, 3. Unless 6 = c = we obtain only the trivial solution. 
Otherwise for f{u) = 

= Kix + n2, T = 2Kit + K3, (f) = —2k\u/ {n — \). (3-14) 
Hence for /ci 7^ (set ki = 1) we obtain the symmetry reduction 

u{x, t) = {t + K3/2)-'/("-') w;(-2), z = {x + K2)/it + K3/2)'/' (3.15) 
where w{z) satisfies 

w" + \zw' + + u;" = 0. (3.16) 

^ n — 1 

It is easily shown that this equation is not of Painleve-type. 

Case 3.5: /(n) = e"" + b-u + c, a 7^ 0. Unless 6 = c = there is only the trivial solution. If 
/(n) = e" then 

,^ = KlX + K2, r = 2Kit + K3, = — 2/Ci. (3-17) 

Hence for 7^ (set ki = 1) we obtain the symmetry reduction 

u{x, t) = w{z) - ln{t + K3/2), z={x + K2)/{t + K3/2)^/2 (3_ig) 
where w{z) satisfies 

w" + i W + 1 + e"' = 0. (3.19) 
It is easily shown that this equation is not of Painlcvc-type. 

Case 3.6: f{u) = ln(n) + bu + c. In this case we obtain only the trivial solution. 

In Table 1 we summarise the infinitesimals obtained using the classical Lie method. For com- 
pleteness, we have included the case when f(u) is linear. 



4 Nonclassical Symmetries 

There are two sorts of nonclassical symmetries, those where the infinitesimal r is non-zero, and 
those where it is zero. In the first case, we can assume without loss of generality that r = 1, while 
in the second case we can assume that ^ = 1. 

In this section, wc solve the determining equations for the nonclassical infinitesimals, and obtain 
the reductions. Solutions are catalogued in §5. 
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Table 1: Infinitesimals for equation (1.1) obtained using the classical Lie method 




1 



u 

(n7^0,l) 
e" 

nln(n) + hu 
arbitrary 



ax + /3t + ^xt + b 
ax + /?t + ^xt + (5 



ax + /3t + ^xt + b 

ax + P 
ax + f3 
ae* + /3 
a 



2at + + K 
2at + + K 



2at + 7*2 + K 

2Q!* + 7 
2at + 7 

7 
/? 



A-i/3x- + 1x2^7 M + 0o(2;,t) 
A - i/3x - (ii + ix2^ 7 u 
+(2a - A)t + ^pxt + (ix2t + |t 



+4>o{x,t) 




X + 2at- \l3x + 


2 1 
2 


+ ^Q{x,t) 




2au 




n — 1 




-2a 




u (^—^ax + 6^ e* 










where cho satisfies 



D,t — 90,xx 



and V'o satisfies -00,* = V'o,xx + V'o 



4.1 Nonclassical Symmetries, r = 1 

The determining equations for the nonclassical symmetries, given r = 1, are 



fl- 


^uu — 


0, 




f2- 


4'uu 


'^ixu + 2ii^ = 0, 




h- 


2(j)xu 


- 2(l)^u + 3^„/ + - (,xx 


+ = 


h- 


4>t - (t 


'^xx + (f>uf - fu4> + 2(A?x - 


- 2^./ = 0. 



These were calculated using the macsyma programme symmgrp.max [16|. 

The result of the "direct search" algorithm on {/i,/2,/3} in the ordering (j) > ^ > f and 
X > t> u \s the equation 

^ufuuuu — 0- 

Hence, unless ^„ = 0, f{u) must be cubic in u. In the case = 0, another "direct search" on the 
system led to the equation 

Cx (,fuu) ^ 2fuu ifuuuu) ~l~ fuufuuufuuuuu ~l~ fuuuu ifuuu) ^ — 

which is the same condition ^3 for f{u) obtained in the classical case. If = = then we obtain 
the condition 

1^2 {fuuu) fuufuuuu^ — 0. 

As in the previous Section, we consider the various possibilities for f{u) separately. 
4.1.1 T = 1, f{u) =u^ + bu'^ + cu + d 

The direct search algorithm yields the equation Cui2^u + 9) = 0. 
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Case 4.1.1i 2^^ + 9 = 0. Here a triangulation of the system yields the infinitesimals 

^ = |i\/2(n + fe/3), (/>= Ku^ + fon^ + cM + d). (4.1) 

and the associated invariant surface condition is 

|i\/2 (u + b/3)u^ +ut- ^{u^ + bu^ + cu + d) = 0. (4.2) 

By solving this we obtain exact solutions (1.1) with f{u) = + bu'^ + cu + d (see §5.1 below). 
Setting ^„ = we obtain the triangulation 

+ + =0, 

6 - se.x + 2ee, = 0, 

•^^xxx '^^^xx ^x (^b 3cj = 0, 
^xx (~3^a;x + 2^^a;) = 0, 

C^x {-b^ + 3c) (3c - 62 - e) (9c - 36^ - 2^^ + 6^^,) = 0, 



Cx (-27d - 2b^ + 96c) = 0. 



Unless the algebraic condition —27d — 26^ + 96c = on the coefficients of f{u) is satisfied, the 
only solution is the trivial one (C = Co a constant, (j) = 0). So assume that 7^ and that 
-27 d - 26^ + 96c = 0. There are thus several cases to consider. Note that if both 6^ — 3c = and 
-27d - 26^ + 96c = then /(n) = {u + 6/3)^. 

Case 4.1.111 = 0, S,xx = 0, f{u) = {u + b/3)^. In this case we have from the triangulation 
equations that Ct + 2S,(,x = and we thus have that 

^ = {x + Ki)/{2t + K2), (p = -{u + b/3)/{2t + K2). (4.3) 



These are equivalent to (3T) and thus also yield the (classical) scaling reduction ([37 
Case 4.1.1111 = 0, —3£,xx + 2^^^ = 0, f{u) = (n + 6/3)^. In this case we have = and 
the Kolchin-Ritt algorithm on the equations {—3(,xx + 2S,(,x = 0, —3S,xxx + 3S,^xx = 0} leads to the 
condition 2S,xi3£,x - C^t) = 0. Thus the non-trivial solution in this case is 

e = -3/(x + Ki), (l) = -3{u + b/3)/{x + Kif. (4.4) 

Hence we obtain the nonclassical symmetry reduction 

u{x, t) = {x + Ki)w{z) — 56, z = \x^ + KiX + 3t (4.5) 

where w{z) satisfies 

w" + u-^ = (4.6) 

which is solvable in terms of elliptic functions (see §5.2 below). 

Case 4.1.11v = 0, 27d + 26^ - 96c = 0, 6^ - 3c / 0. In this case we have 

= 0, 9c - 36^ - 2f + 6^^ = 

leading to the solution 

C = 3//tan(/ix + Ki) = — /i(3n + 6) sec2(/ux + Ki) (4.7) 
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where /u^ = — 3c), and so we obtain the nonclassical symmetry reduction 

u{x,t) = fj,cot{iJ,x + Ki)w{z) — ^b, z = sin(^x + Ki) exp(— 3/i^t) (4-8) 
where w{z) satisfies 

z^w" - 2zw' + 2w + w^ = (4.9) 

which is solvable in terms of elliptic functions (see §5.2 below). 

We remark that the roots of the cubic + bu^ + cu + (9c — 26^)6/27 = are —^b, —^b + V2 n 
and — \/2/i, where fi"^ = g(6^ — 3c), i.e. they are collinear and the distances from the outer roots 
to the central one are equal. The condition that 27d + 2b^ — 96c = is the necessary and sufficient 
condition on f{u) = + bv?' + cu + d for the existence of nonclassical symmetry reductions of (1.1) 
to ODES solvable in terms of elliptic functions. Furthermore, Nucci and Clarkson |71] show that 



elliptic functions solutions exist for the Fitzhugh-Nagumo equation (|l.9D if either a = — 1, a = ^ 
or a = 2, which are precisely those the three cases when the distances from the outer roots to the 
central one are equal. 

4.1.2 T = 1, f{u) = au + b 

A preliminary triangulation yields the condition ^„ = 0. Reducing the equations with respect to 
= we obtain = 0. Inserting cf) = A{x, t)u + B{x, t) into the equations, we obtain the system 

in 

gi : 2A, + 6 - + 2^?. = 0, 

52 : At- A^^ + 2 {A -a) = 0, 

93 ■■ Bt- B^^ + 2{B-b)^^-aB + bA = 0. 

Note that we no longer have an over-determined system. For an arbitrary function v{x,t) let 
^{v; C) = ~ Vxx + 2v(,x, and set $(,^; ,^) = r/. Let A > ^. Then cross-differentiating gi and g2 and 
reducing, we obtain 

A^xx(.A-a) = ^rj;0- 

Case 4.1.21 r = 1, ^^x / 0, /(n) = au + b. Pseudo-reducing gi, g2 with respect to ^xx{A — a) = 
j^{r];^), we obtain two equations ki and ^2 in C whose highest derivative terms in a total degree 
ordering are ^xxxxxx and ^xxxxx- (A total degree ordering is determined first by the degree of 
the terms and then by a lexicographic ordering.) Reducing ki with respect to k2 in this ordering 
yields a condition k^ for ^ whose highest derivative term is S,xxxxt- Now {^2,^3} forms a DGB for 
I{{ki, k2}) with 5({/c2, ks}) = {± (^^x)" : n £ N}. (We have stipulated that S,xx 7^ in this case.) 
Our next step is to apply the Initial Data algorithm |Q] , from which we have that a formal solution 
to the equations for ^ depends on four arbitrary functions and one arbitrary constant. In fact we 
can assign ^(0, t), ^x{0,t), Cxx{0,t), CxxxiO,t) and ^xxxxi^^O) and the other coefficients for a formal 
expansion for ^ will be determined. 

Assuming that ^a:a:(0, i) is a constant, and that the other arbitrary functions and constant are 
zero, one can derive a solution to the equations for ^ and A (i.e. gi and (72): 

^ = T^' ^ = a + -^ (4.10) 

3c — x-^ 3c — x-^ 

where c is a constant, simply by summing the formal Taylor series obtained from k2 and k-^. Solving 
53 for B with ^ and A as given in ( [4.10D yields 

B — b+ "^^^^ 

a(3 — cx"^) 
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Therefore the infinitesimals in this case are 

, , 3x(au + b) 

3' <P={au + b) + —\ ^ (4.11) 

6c — x'^ a{6c — x'^) 

and so we obtain the symmetry reduction 

u{x,t) = xw{z)e'^^ -b/a, z = \x^ + t + c/x (4.12) 

where w{z) satisfies 

w" = 0. 

Hence we obtain the exact solution of (1.1) with f{u) = au + b given by 

u{x, t) = \a{\x^ + xt + c) + /?} e"* - b/a. (4.13) 

Case 4.1.211 r = 1, ^^x = 0, f{u) = au + b. In this case we obtain a classical reduction and so 
we omit details. 

Nonclassical symmetry reductions for the linear equation 

ut = Uxx + au + b (4-14) 

where a and b are arbitrary constants, obtained using both Painleve and symmetry analysis of the 
determining equations {gi, 92,93} can be found in [98, pS| ]. 

4.1.3 Nonclassical, r = 1, f{u) neither cubic nor linear 

In this case we have 

'^fuu {fuuuu} + fuufuuufuuuuu + fuuuu ifuuu) — 



or 



2 {fuuu) fuu fuuuu — 0- 



The various possibilities for f{u) have already been listed in §3. For all these possibilities, we obtain 
the same symmetries as for the classical case, with the classical ^ /r yielding the nonclassical ^ and 
the classical cj) /t yielding the nonclassical (p. However, for some of the cases, this is harder to prove 
than for the classical determining equations. 

We give the proof for the case f = u^ + cu + d where n 7^ 0, 1, 2 or 3, to show the ideas involved. 

We have in this case that = 0. The Direct Search algorithm yields that d = else one has 
only the trivial solution, and that (f) = —2u^x/{n — 1). 

Reducing the determining equations with respect to and + 2 /(n — 1) uS,x yields 

hi : ^ixxx - 2(n - \)iixx + c{n - Xfi^ = 0, 
^2-- in- \)it - (n + 3)^^^ + 2(n - 1)^^^. = 0. 

Take the ordering t > x. Then cross-differentiating hi and /i2, throwing away powers of c, and 
(n — 1) that factor out, and iterating, we obtain successively equations with the following "vital 
statistics" : 

/13: HDT(/i3)^P('^^) = exx, Sep(/i3) = 24<.. + n\n - 
hi-. HDT(/i4)Hp('^*) = Cxx, 

Sep(/i4) = -12(n + 1)^^ + (5n2 - 2n^ - 5)^^ + 12c(n2 - 120n + 108), 
/15: HDT(/i5)^P('^^) = e^, Sep(/i5)=p(C.,0 
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where p is a polynomial. Using the Implicit Function Theorem, we obtain for some function R, the 
ansatz = R{0- Note that if any of the separants are zero, leading to invalid reductions, we still 
obtain the same ansatz. Substituting the ansatz into hi, /12, one obtains 



h 



- R{() = 0, 
(n - l)Ct + 2(n - 1)CRR^ - (n + 3)RR^ = 0, 

(n - + 2(n - 1)CRR^ - ARR] - 2{n - l)R^ - n{n - 2)cR = 0. 



Cross-differentiating ki and /c2 > R) and reducing yields 

R^ [{n + 3)% - 2{n - 1)] = 0. 



(4.15) 



Solving for R using the second factor of (|4.15D , and substituting the result into k^, yields a poly- 
nomial for ^, that is, ^ is a constant. 

Thus, since the only factors we threw away were c, and (n — 1), it must be that for a nontrivial 
solution to exist, c = 0. 

Setting c = in /ii and /i2 and performing the Kolchin-Ritt algorithm on these two equations 
yields the condition 



n 



The system {^^ 



(n - i)e 

0, + 2^(^x = 0} has for its solution 

X + K2 , 



mx - i2Cx 



(4.16) 



2u 



2t + Ki' (n-l)(2t + Ki) 

which is equivalent to the classical reduction of Case 3.3. 

Putting the other factor of (4.16), namely (n^ — l)S,(,x — ^2^xx, in addition to hi and /12 (with 
c = 0) as input to the Kolchin-Ritt algorithm, we obtain 

- l)(n + 1) [12e.. + (n - l)(n - 5)^] = 0, 
e^(n-l)3(n-3)(n- 5)^ = 0. 



Since n 7^ 1 or 3 in this case, there are no other nontrivial solutions. 



4.2 Nonclassical Symmetries, r = 0, ^ = 1 

We come now to those nonclassical symmetries where r = 0, and" without loss of generality we can 
take ^ = 1. In this case we have the system consisting of the equation (1.1), the invariant surface 
condition (1.4) and one additional determining equation for <f>: 

/i : ut- Uxx - f{u) = 0, 
/2 : Ux - <l)ix,t,u{x,t)) = 0, 

/s : (j)t- (t>xx - 4>'^(t>uu - 2(f)(t)xu + f<t>u - fu(t> = 0. 

The last equation was calculated using the macsyma package symmgrp.max Cross-differen- 
tiating /i and /2 produces exactly /s so we obtain no new conditions. So we try three ansatze for 
0, namely 

Ansatz 4.2.1: (j) = P^u^ + Piu + po with f{u) = —v? — hv? — cu — d. 
Ansatz 4.2.2: = gi(x, t)u + qQ{x, t). 
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Ansatz 4.2.3: cf) = ri{t,u)x + rQ{t,u). 

In the second and third ansatze, f{u) is to be determined. We consider these three ansatze in turn. 
Note that in the equation /s, u is to be regarded as an independent variable. 

Ansatz 4.2.1. Reducing using the ansatz, and setting coefficients of powers of u to zero, we 
obtain from the coefficient of that p2 {2p2 — 1) = 0. Setting p2 = we have from the coefficients 
oi v} , u and 1 that: 



91 

92 

93 



2pl - V2bpi - V2po + c = 0, 
2poPi - V2bpo + d = 0, 
V2pl - cpo + dpi = 0. 



The coefficient of is identically zero. For generic b, c and d there are three possible solutions of 
{91,92,93} given by 

pi = -^V2{mi + nij), po = ^V2mimj, l<i<j<3, (4.17) 
where rrii and mj are any two of the three distinct roots of the cubic 

+ bu'^ + cu + d = 0. (4.18) 

Hence there are three choices for (f> 

(u — mi){u — rrij), 1 < * < J < 3. (4.19) 
Substituting this into /2, solving for u and requiring that u also satisfies /i yields the solution 

ut = Uxx — — bu^ — cu — d = Uxx — {u — mi){u — m2){u — ms), (4.20) 

given by 

t) = ' ' ^ '\ \ l<^<3<3, 4.21 

where 

0fc(x, t) = exp |i\/2mfcX — (jrii + m,2 + — |mfc^ t| , 1 < A; < 3. 

We remark that the solution ( [4.21D is special case of the travelling wave solution and of the expo- 
nential solution ( ^.13| ) given in §5.1.1 below. 

If d = (9c — 26^)6/27, then there is the solution of {91,92,93} given by pi = |\/2 6 and pQ = 
^\/2{c— §6^). Thus, a similar procedure to the generic case above yields exact solutions of 



ut = Uxx -u^ - bv? - cu - (9c - 26^)6/27, 



given by 



(c - i62) 1/2 tan (i^/2 (c - |62) 1/2^; + _ 1 5^ jf ^ ^ 1 6^ 
ift, if c = ift^ 



X + K2 

where ki and At2 are arbitrary constants. We remark that d = (9c — 26^)6/27 is the same relation 
as that we obtained for (1.1) with f{u) cubic and r = 1 which yielded special reductions in terms 
of elliptic solutions. 



Symmetry Reductions of a Nonlinear Heat Equation 



18 



If pi = then po = ^ \/2 c provided that d = be. Hence we obtain the exact solution of 

ut = Uxx — — hu^ — cu — be, 

given by 

u{x,t) = ^/c tan \/c (^\/2 X — 6t) , 

which is a special case of the travelling wave reduction. 

If p2 = 0, we then have that pi = po = implying that (p = 0. Thus solving the equations fi 
and /2 yields 

u(x, t) = g{t), where ^ + ag^ + bg^ + eg + d = 0. 

dt 

Ansatz 4.2.2. In this ansatz, we take 

(j) = qi{x,t)u + qo{x,t) (4.22) 

Reducing using the ansatz, it can be shown that there are three cases to consider. 

Case 4.2.21 / = au^ + 6ti + c, a / 0. In this case we have cj) = and solving the equations /i and 

/2 yields 

, , , , dg 2 7 
u[x,t)=g{t), +bg + e. 

Case 4.2.2ii f = au + b. In this case, it is easily shown that </> = (u + b/a)qi, where qi satisfies 

qi,t - qi,xx - 2gi,xgi = 0. 
Combining this with the equations /i and /2, we recover the well-known Cole-Hopf transformation 



22, 51 1 relating a solution of Burgers' equation to a solution of the linear heat equation. 



Case 4.2.2iii f{u) = {u + k) ln{u + k) — {u + k) + du + e. In this case, it is easily shown that if 
e = dk then = (n + k)qi{x, t) where qi satisfies 

qi,xx + qi- qi,t + '^qiqi,x = 0, 

whilst if e ^ dk then = 0. The situation is analogous to that in Case 4.2.2ii above. 

It appears that Backlund transformations arise in considering nonclassical symmetries in the 
r = case. Further interesting examples have been obtained by Nucci |70|. 

Ansatz 4.2.3 (p = ri{t,u)x + rQ{t,u). Reducing using the ansatz, it can be shown that there 
are three cases to consider: (i), /(n) = au^ + ftu + c, a 7^ 0, (ii), f{u) = au + b, a ^ 0, 01 (iii), 
f{u) = u ln{u) — u + eu + d. 

Case 4.2.3i f{u) = au^ + bu + e, a 0. In this case, (j) = 0. 

Case 4.2.3ii f{u) = au -\- b, a ^ 0. In this case we obtain the exact solution of 

Ut = Uxx + au + b, (4.23) 

given by 

, Ki { {x + Kif\ b + Kse"^ 

u[x,t) = — 77rexp<at — — ^> , (4.24) 

(i + K2)V2 ^1 4(t + K2) J a ' ^ ^ 

where ki,K2,K3 and K4 are arbitrary constants. 
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Table 2: Infinitesimals for equation (1.1) obtained using the nonclassical method 



/(«) 



+ bu^ + cu + d 
—u^ — bu^ — cu — d 

+ hv? + CU + 6(9c - 26^) /27 
-u^ - bu^ -cu- 6(9c - 262)/27 

au + b 

— {u — mi){u — m2){u — ma) 

— (u — mi)('u — m2){u — 7713) 

— (u — mi){u — m2)('u — 7713) 
-u^ - \m? -cu- 6(9c - 252)/27 

— — — cu — be 
{u + A;) ln(?i + k) + c{u + A;) 



iV2(u + 6/3) 
V2(u + 6/3) 
3 

X + Kl 

3 



3;utan(/xx + K3) 
3/xtanh(/xx + K4) 
3x^ 

3^5 — 

1 
1 
1 
1 
1 
1 



'^{u^ + bu^ + cu + d) 
-^{u^ + bu^ + cu + d) 
3(n + V3) 

3(n + 6/3) 

(X + Ki)2 

— )u(3u + b) sec^(/xx + K3) 

/i(3u + b) sech^(/xx + /«4) 

3x(au + b) 

au + b^ — — ^T- 

a{6K^ — x^) 

^\/2 (tx — mi){u — 777.2) 
^y/2{u — 7772)(w — 7773) 

^\/2 (u — 7773)(u — 777l) 



hV2 



[u + b/^f + (c - \b^] 



\^/2{u^ + c) 



-\{x + KQ){au + 6) 

(x + K8){au + b + Kge"*) 
2(t + Kio) 



where ki, K2, ■ • ■ , /^lo are arbitrary constants. 

Case 4.2.3iii f{u) = {u + k) ln(7i + k) + c{u + k). In this case we obtain the exact solution of 

ut = Uxx + {u + k) \n{u + k)+ c{u + A;), (4.25) 

given by 



u{x,t) = exp |K2e* — ^(x + ki)^ — c + ^| — A; 
where ki and K2 are arbitrary constants. 

5 Exact solutions of Ut = Uxx + + hv? -\- cu-\- d 

In this section we obtain exact solutions of the pdes 

Ut = Uxx + u^ + bu^ + CU + d 

and 

■"t = Uxx — + bu^ + CU + d) 



(4.26) 



(5.1) 



(5.2) 



using the infinitesimals derived in §4.1.1 above. There are two types of solutions: those expressible 
in terms of exponentials arising from Case 4.1.1i (§5.1) and those expressible in terms of elliptic 
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functions arising from Cases 4.1.1iii and 4.1.1iv (§5.2). It is interesting to note that the form of the 
exponential-type solution is closely related to the roots of the cubic 

+ bu^ + cu + d = (5.3) 



arising in the pde we are considering. Furthermore the conditions when (5.1,5^) possess solutions 



expressible in terms of elliptic functions are determined by the roots of ( ^.3| ) . Recall that nonclassical 
symmetry reductions of (5.1), (5.2) are obtainable in the case when the roots of f{u) are arbitrary 
(Case 4.1.1i), and when the roots of f{u) are collinear and the distances from the outer roots 
to the middle root are equal (Case 4.1.1iv). These reductions lead to exponential and elliptic 
type solutions respectively. The second case is determined by the equivalent algebraic condition 
27 d + 2b^ - 96c = on the coefficients of f{u) = + hv? + cu + d. We note that a solution of (O) 



can be obtained from one of (5J) by letting x — > itix and t —t, and vice- versa. 



5.1 Exponential-type solutions 

In Case 4.1.1i above, we obtained the infinitesimals for ( ^.1[ ) given by 

C = |i\/2(u-F 6/3), (p = ^{u^ + bu"^ + cu + d) (5.4) 
with associated invariant surface condition 

|i\/2 {u + b/3)ux +ut- ^{u^ + bu^ + cu + d) = 0. (5.5) 



Following Nucci and Clarkson ||7l|], this is solvable as follow. Using (^^) to eliminate Ut in ( p.lj) 
yields 

u^x + |iV2 (u + b/3)ux -^{u^ + bu^ + cu + d) = Q 

which can be linearized by the transformation u = — \/2i*&2:/*^ yielding 

2^/2^xxx + 2ib^xx-V2c^x-''^d<^ = ^■ (5.6) 

Suppose that the roots of the cubic (|5.3[) are mi, m2 and ms, then the roots of 

2V2p^ + 2\bp^ - ^/2cp-id = {) (5.7) 

are ^\/2imi, ^\/2im,2 and ^\/2im^. 

The analogous infinitesimals for (|5.2[) are 



C = |\/2(u-F6/3), (I) = -^u^ + bu'^ + cu + d) (5.8) 
with associated invariant surface condition 

^V2{u + b/3)ux + ut + l{u^ + bu^ + cu + d) = 0. (5.9) 
Using this to eliminate Ut in (|5.2| ) and making the linearizing transformation u = \/2^'a;/^ yields 

2V2^xxx + 2b^xx + V2c^x + d^ = 0. (5.10) 

The roots of 

2V2q^ + 2bq'^ + V2 cq + d = (5.11) 

are ^V2mi, ^V2 m2 and |\/2m3, where mi, m2 and ms are the roots of the cubic ( p!3|) . 

In the following four subsections we consider the cases when mi, m2 and ms are real and distinct, 
complex, two are equal and all three are equal, respectively. These exponential-type solutions of 



|5.2| ) are real for real x and t whereas those of (5.1) are complex for real x and t. Consequently 



we only give details for ( |5.2| ) and leave it for the reader to derive the analogous solutions for (5.1) 
using the transformation x ibix, t — > —t. 
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5.1.1 Distinct Real Roots 

Solving ( |5.10| ) in the case when mi, 1712 and are distinct and real yields 

^{x, t) = ^i{t) exp mix^ + /U2(t) exp {^V2m2x) + l^3{t) exp msx^ (5-12) 

where fJ-i{t), fJ-2{t) and fis{t) are to be determined. By requiring that u = \/2$a;/$ also satisfies 
(|5.9|), or equivalently (^), it is easily shown that 



(t) = Kj exp |m,j (jni + m2 + ms — ^rrij) t| , j = 1, 2, 3, 



where ki, K2 and K3 are arbitrary constants (note that mi + m2 + ms = —6). Hence we have 
obtained the following exact solution of ( ^.21 ) given by 



_ Kimi^i(x, t) + K2m2^2(a^, t) + K3m3^3(x, t) 

" Kl^l(x,t) + K2^'2(x,t)+K3^3(x,t) 

where 

t) = exp ||\/2 mjX — rrij {mi +m2 + m^ — fmj^ t| , j = 1, 2, 3. 

Setting mi = a, m2 = 1 and m^ = in ( 5.13 ), which we may without loss of generality for 
the general cubic with real and distinct roots, yields the following exact solution of the Fitzhugh- 
Nagumo equation (|1.9| ) for a ^ and a 7^ 1 



ani exp 1 i {±^2 ax + a^t 


} + «;2exp|i(±^/2x + t^ 


} 


Ki exp 1 i 


'±V2ax + aH 


} + K2exp{i (±^/2x + t^ 


1 + K3 exp {at) 



^(^' ^) = rw m rw M • (5-14) 



This solution was obtained by Vorob'ev (who calls the associated symmetry a "partial symmetry 
of the first type"), Kawahara and Tanaka |^5| (using Hirota's bi- linear method |50| ]) and Hereman 



||45|| (using the truncated Painleve expansion method ||67| , 100 , 101]). Plots of ( 5.14 ), with ki = 
1-^2 = = Ij for (i), a = —0.5, (ii), a = 0.4, (iii), a = 0.7, and (iv), a = 1.5 are given in Figure 1; 
these plots were drawn using maple. 



5.1.2 Complex Roots 

Suppose in this case that the roots of ( p.3| ) are m (real) and a ib i/3. Then setting mi = m, 
777-2 = a + i/3 and 7773 = a — i/3 in ( ^.13|) we obtain the exact solution of 

ut = Uxx - + (777 + 2a)7i^ - (a^ + + 2a777)7t + (a^ + (3'^)m (5.15) 

given by 



Ki7n,exp l^-v/2/ia; + 


+ 


+ K2 (a^ + /32)i/2 (^^y/2px - nl3t 


+ O0 + 60) 


Ki exp ij,x + (^^/U^ 




t| + K2 cos 


[^^V2 px - fif3t + So) 





(5.16) 

respectively, where ki, K2 and Sq are arbitrary constants, = m — a and = tan ^(/?/a). 
Setting 777 = a = and /3 = 1 in (|5.16| ) yields the exact solution of 

Ut = Uxx - u{u^ + 1) (5-17) 
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given by 



u{x, t) = 



K2 sin (^^\/2x 




Ki exp (^^t 


+ K2 cos (^^V2x^ 



Setting m = and a = /3 = 1 in ( 5.16| ) yields the exact solution of 

ut = Uxx — u{u^ — 2u + 2) 

given by 

u{x, t) 



2; these plots were drawn using maple. 
5.1.3 Two Equal Roots 



K2 


COS ^|\/2 X — t 


+ sin(i^/2x-t) 




Ki exp 

K2 = 1, 


'\'sf2x + 2i) + K2COS 
and (|5.20|), with k\ = 


[^V2x-t) 
1 and H2 = 



(5.18) 



(5.19) 



(5.20) 



jure 



Suppose that the cubic ( ^ ) has a single root mi and a double root m2. Then the solution of ( ^.10 ) 
is given by 

^'(x,t) = ^i(t)exp (^i\/2mix) + [n2{t) + x/X3(i)] exp (i\/2 m2x) (5.21) 

where /ii(i), /^2(i) and //3(t) are to be determined. A similar procedure to that used in the previous 
subsection yields the exact solution of 



Ut = Uxx - {u- mi){u - m2f 



given by 



m(x, t) 



niiKi exp 1^ 


'V2f3x + p^i^ 


}+K2 


1712 (x - V2Pt 




Ki exp 


{i [V2px + I3H) 


}+K2 (^X-V2pt) 



(5.22) 



(5.23) 



where (3 = m2 — mi and ki, K2 and K3 are arbitrary constants. 
In particular, an exact solution of 



Ut = Ux 



u^{u + h) 



IS 



u{x, t) 



bni exp 1 i 


'y/2bx + b'^i 


} + ^/2k2 


Ki exp |i (y2bx + bH 


} +K2 


'x - y/2bt^ 



(5.24) 



(5.25) 



Plots of this solution for b = ±1, ki = 1 and ^2 = ^ are given in Figure 3; these plots were drawn 

using MAPLE. 

5.1.4 Three Equal Roots 



In the case when the cubic ( p.Sp has three equal roots, a similar procedure to that used in the 
previous three subsections yields the exact solution of 



Ut = Uxx — U^ 



(5.26) 
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Table 3: Exponential- type solutions of equation (1.1) 

u{x, t) 



u{x, t) 

aKi exp I i (^±\/2ax + aH^ | + ^2 exp | ^ [±^/2 a; + t j | 

Hi exp I i (^±\/2 aj; + a^tj | + K2 exp | ^ ^ib-v/2 x + | + K3 exp 

K,o sm ^i^/2 



ii(l — ti)(u — a) 

-u{u^ + 1) 

-n(ti2 -2^ + 2) 

-u^{u + 6) 
— n'^ 



a; 



it) 



K2 sin 


[\V2x] 


) 


Ki exp {^t^ + ^2 cos 1 


[\V2x) 



K2 


cos 


'^V2x-t 


) + sin (^i\/2x - t'j 




Ki exp (^^\/2x + 2t^ 


+ K2 COS (^^\/2 X — 



6ki exp |i 


'V2bx + bH^ 


} + \/2k2 


Ki exp 1 i 


'y/2bx + hH 


} +K2 (x-^/2 6t) 



+ KiX + Qt + K2 



given by 



u{x, t) 



V2{2x + Ki) 

+ KlX + Qt + K2 



(5.27) 



A plot of this solution for ki = K2 = is given in Figure 4; this plot was drawn using maple. 

In Table 3, we list the nonclassical exponential- type solutions of equation (1.1) derived in this 
section for various canonical choices of f{u). 

5.2 Elliptic Function Solutions 

In Case 4.1.1iii above we derived the nonclassical symmetry reduction of 

,3 



given by 
where w{z) satisfies 
The solution of this equation is 



Ut = Uxx + U 



u{x,t) = {x + Ki)w{z), Z = ijX + KiX + 3t 



w{z) = i\/2 sd {z;\V2 

where sd(z; k) is the Jacobi elliptic function satisfying 

2 



^) =l + (2A:2-l)r?2 + fc2(fc2_i)^4^ 

Hence we obtain the exact solution of (|5.2^ ) given by 

, t) = i^/2 (x + Ki) sd [\x^ + Kix + 3t; i V2 



n X, 



(5.28) 

(5.29) 
(5.30) 
(5.31) 

(5.32) 
(5.33) 



A plot of this solution for ki = is given in Figure 5; this plot was drawn using mathematica. 
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The analogous nonclassical symmetry reduction of 



Ut = 



is given by 
where w{z) satisfies 
The solution of this equation is 



u{x, t) = {x + Ki)w{z)^ Z = + KiX + 3t 



w" - = 0. 



w{z) = \/2 ds (z;i\/2 
where ds(2;; A;) is the Jacobi elliptic function satisfying 

(^^y = A;2(fc2-l) + (2A;2-l)r?2 + r?4. 

Hence we obtain the exact solution of ( |5.34| ) given by 

u{x, t) = V2{x + Ki) ds [\x^ + Kix + 3t; \V2j . 

In Case 4.1.1iv above we derived the nonclassical symmetry reduction of 

6(9c - 262) 



Ut = + u +hu + cu + 
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(5.34) 

(5.35) 
(5.36) 
(5.37) 

(5.38) 
(5.39) 
(5.40) 



given by 



u[x 



, t) = Ki/i sin(^x + K2) exp ^— S/i^t^ w{z) — ■^h, z = ki cos(;UX + K2) exp (^—S^'^t^ (5-41) 



where /i'^ = ^{b'^ — 3c), ki and K2 are arbitrary constants and w{z) satisfies (5.30). Hence for > 3c 
we obtain the exact solution of ( |5.40| ) given by 

ti(x,t) = i\/2Ki/isin(^a; + K2) exp (—3/^2^1 sd Kicos{fix + K2)exp (—Sfi^t) -,^^2 — hb (5.42) 



where jj"^ = ^(^^ — 3c). Two plots of this solution for ki = 3, K2 = and /x = ^ are given in Figure 
6; these plots were drawn using mathematica. For 6^ < 3^ we obtain the exact solution of ( p.40| ) 
given by 



u[x, 



, t) = |\/2 KiA sinh(Ax + K2) exp |^3A2t^ sd ki cosh(A3; + K2) exp ^SA^t^ ; i\/2 — ^b (5.43) 



where A^ = i(3c — 6^). A plot of this solution for ki = 3, k;2 = and A = ^ is given in Figure 7; 
this plot was drawn using mathematica. 

The analogous nonclassical symmetry reduction of 



, , 2 6(9c-262 
Ut = Uxx — u — bu — cu — 
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(5.44) 



is given by 

u{x, t) = Ki/x sinh(/ij; + K2) exp (Sfi^t^ w{z) — ^b, z = cosh.{fix + K2) exp (^3^'^t^ (5.45) 
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Table 4: Elliptic Function solutions of equation (1.1) 



u(x, t) 



u{u^ - 2) 
n(ti2 + 2) 
-n(ti2 - 2) 
-u(u'^ + 2) 



i\/2 (x + Ki) sd (ix^ + Kix + 3t; i\/2 
V2{x + Ki) ds (ix^ + Kix + 3t; ^V2j 
i\/2 Ki sin(x + K2) exp (— 3t) sd ki cos{x + K2) exp (— 3t) ; ^\/2 
^\/2 Ki sinh(x + K2) exp (3t) sd ki cosh(x + K2) exp {3t) ; ^\/2 
\/2 Ki sinh(x + K2) exp (3t) ds ki cosh(a; + K2) exp {3t) ; ^\/2 
\/2 Ki sin(x + K2) exp (— 3t) ds ki cos(x + ^2) exp (— 3t) ; ^\/2 



where //^ = — 3c), ki and K2 are arbitrary constants and w{z) satisfies (5.36). Hence for 6^ > 3c 



U (5.46) 



we obtain the exact solution of ( ^.44 ) given by 

u{x,t) = \/2 Ki/i sinh(/ia; + K2) exp (^3/i^tj ds ki cosh(^a; + K2) exp (S/i^t^ ; 
and for 6^ < 3c, obtain the exact solution 

u(x, t) = \/2 KiAsin(Ax + K2) exp 3A^tj ds ki cos(Ax + K2) exp ^— 3A^t^ ; ^\/2 — |6 (5.47) 
where X"^ = l{3c - b"^) . 

In Table 4, we list the nonclassical elliptic function solutions of equation (1.1) derived in this 
section for various canonical choices of f{u). 



6 Discussion 

In this paper we have demonstrated that the method of dgbs has enabled all symmetries of the 
nonlinear heat equation (1.1) to be found; the list it provides of possible analytic functions f{u) for 
which symmetries may exist is definitive. The use of dgbs has made the analysis of overdetermined 
systems of pdes more tractable and whilst the diffgrob2 package needs to be used interactively at 
present, nevertheless it is effective in solving the overdetermined systems of determining equations 
for classical and nonclassical symmetries of (1.1). 



It is not clear how the direct method, developed by Clarkson and Kruskal |20| for finding 
symmetry reductions of PDES may be applied to equations such as (1.1) which contain arbitrary 
functions. Nucci and Clarkson [^l| (see also |21]) demonstrated that the symmetry reduction ( ^.14 ) 
of the Fitzhugh-Nagumo equation ( |1.9|) , is obtainable using the nonclassical method due to Bluman 
and Cole |^], though not using the direct method due to Clarkson and Kruskal Hence the 

nonclassical method is more general than the direct method, at least as it was originally formulated. 



Olver 1 74] has recently shown that the direct method is equivalent to the nonclassical method when 
the infinitesimals and r are independent of the dependent variable u (see also [19, 81, 104, 105] ). 
In this case the associated vector field is of the form 



w = ^(x, t)dx + t(x, t)dt + (/)(x, t, u)di 
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and generates a group of ^^fibre-preserving transformations^^ , meaning that the transformations in 
X and t do not depend upon u. A recent extension of the direct method by Estevez [^] does yield 
the symmetry reduction solution ( |5.14| ) of the Fitzhugh-Nagumo equation ( |1.9D . These results 
pose the following important open question: "for which PDES does the nonclassical method yield 
more symmetry reductions than the direct method?" Furthermore, it remains an open question 
to determine a priori which pdes possess symmetry reductions that are not obtainable using the 
classical Lie group approach. 

Fushchich and Serov (see also ^3|) claimed that 

ut + Uxx = fiu) (6.1) 

which is equivalent to (1.1), has a nonclassical symmetry, which they refer to as a conditional 
symmetry with respect to the vector field (1.3), if and only if (1.1) is equivalent to the special form 

ut + Uxx = + bu + c (6-2) 

where a, b and c are arbitrary constants, which is equivalent to (5.1) and (5.2) and that the 
associated infinitesimals are 

(, = lV2au, T = l, (j)= ^{au^ + bu + c). (6.3) 

Subsequently they consider the four special cases of (|6.2[) given by 



ut + Uxx = au{v?-l), (6.4) 

Ut + Uxx = au{u^ + 1), (6.5) 

Ut + Uxx = a{u'^ — 3u + 2), (6.6) 

Ut + Uxx = au^ (6.7) 

and whereas they integrate the invariant surface conditions, they do not write down the associated 
exact solutions of (|6.1|) . 

The results we have obtained in this paper show that there are more nonclassical symmetry of 



(6.1) other than those generated by the infinitesimals (|6.3| ). In particular they do not obtain the 
nonclassical symmetry reductions corresponding to those obtained in Cases 4.1.1iii,iv and 4.2.3ii,iii 
above. Further even for the symmetry reductions associated with (|6.3|), the four "canonical" equa- 



tions (|(L4,p3,pl^j6.7|) do not include the case when the roots of the cubic au + bu + c are real 



distinct and unequally spaced. 
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Figure Captions 

Figure 1: Exponential solution ( |5.14| ) of ut = Uxx + ^^(1 — u){u — a) where (i) a = —0.5, (ii) 
a = 0.4, (iii) a = 0.7, (iv) a = 1.5. 

Figure 2a: Exponential solution ( ^.18| ) of ut = u^x — u{v?' + 1) where ki = K2 = 1. 
Figure 2b: Exponential solution ( ^.20D of ut = Uxx — u{v? — 2n + 2) where ki = K2 = 1. 
Figure 3: Exponential solution ( ^.25| ) of ut = Uxx — u{v? + h) where (i) h = —1 and (ii) 6 = 1. 
Figure 4: Rational solution ( |5.27| ) of ut = Uxx + 'u? where ki = K2 = 0. 
Figure 5: Elliptic solution ( |5.33| ) of ut = Uxx + where ki = 0. 

Figure 6: Elliptic solution (|5.42| ) of ut = Uxx + v? + cu + 5(9c - 26^) /27 where = ^{iP' - 3c) = 
Ki = 3 and K2 = 0. 

Figure 7: Elliptic solution (|5.43| ) of ut = Uxx + + cu + 6(9c - where = i(3c - 6^) = 

^\/2, Ki = 3 and K2 = 0. 



